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The inelastic scattering of electrons which carry current through a single-molecule junction is 
modeled by a quantum dot, coupled to electron reservoirs via two leads. When the electron is on 
the dot, it is coupled to a single harmonic oscillator of frequency lvq- At zero temperature, the 
resonance peak in the linear-response conductance always narrows down due to the coupling with 
the vibrational mode. However, this narrowing down is given by the Franck-Condon factor only for 
narrow resonances. Contrary to some claims in the literature, the linear-response conductance does 
not exhibit any side-bands at zero temperature. Small side-bands, of order exp[—f3huJo], do arise 
at finite temperatures. The single-particle density of states exhibits discontinuities and logarithmic 
singularities at the frequencies corresponding to the opening of the inelastic channels, due to the 
imaginary and real parts of the self-energy. The same singularities also generate discontinuities and 
logarithmic divergences in the differential conductance at and around the inelastic thresholds. These 
discontinuities usually involve upwards steps, but these steps become negative within a rather narrow 
range of the elastic transparency of the junction. This range shrinks further as the the excitation 
energy exceeds the bare resonance width. 
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I. INTRODUCTION 

Single-molecule junctions based on direct bonding of a 
small molecule between two metallic electrodes seem by 
now rather established experimentallyFii^i^ii^i^i^iii^i^ The 
electronic transport through such a molecular bridge is 
attracting a great deal of interest, including the inven- 
tion of ingenious experimental realizations for it (see, for 
example, the recent Refs. [l^)- Besides the possible tech- 
nological advantages of "molecular electronics",^^ there 
are many issues that make this problem of great interest 
from both the basic science and the application points 
of view. The possibility of directly addressing a single 
microscopic quantum system with an ordinary measure- 
ment apparatus should shed light on fundamental quan- 
tum measurement questions. 

Electrons passing through the small molecule may 
change its quantum state (electronic, vibrational, and in 
certain cases also rotational, and even the conformatior>i2, 
of the molecule). These may require a finite energy 
transfer from the transport electron. Thus, the dynam- 
ics of the molecule may create interesting structures in 
the I-V characteristics ! ^'^'^'^ These rich characteristics, re- 
sembling the one observed in inelastic electron tunneling 
spectroscopy (lETS),-^ depend on important experimen- 
tal details such as the equilibration time of the vibra- 
tions compared to the typical time between consecutive 
electrons passing through the junction, or whether the 
electrons can pump the molecule into higher vibrational 
states. Such measurements provide a handle on study- 



ing molecular properties and their modifications by the 
binding to the electrodes. In some cases they may also 
help to identify the molecule which has been bound in 
the bridge. 

The configurational modification of the molecule by 
the tunneling electron is usually described by a linear 
coupling of the electron with e.g. the vibrational modes, 
while the oscillating location of the whole molecule is 
modeled by the dependence of the tunneling matrix 
elements to the leads on the vibrational degrees of 
freedomJ^ As is well-known, one may eliminate the lin- 
ear electron-phonon interaction by a canonical transfor- 
mation which dresses the tunneling matrix elements by 
the phonon cloud (the Holstein polaron).'^ The result- 
ing matrix elements contain the Franck-Condon factors. 
These tend to block the conductance at off-resonant situ- 
ations (the Franck-Condon blockade) However, the 
top of the resonance conductance is not reduced by these 
factors. Indeed, we find that the coupling to the vibra- 
tional mode causes a narrowing of the resonance. This 
narrowing is described by the 'usual' Franck-Condon 
blockade only in the limit of very narrow 'bare' reso- 
nances, namely very long dwell times of the electrons on 
the resonances. We find that reducing this dwell time 
weakens the Franck-Condon blocking. 

Transport through small molecules offers new means 
of studying tunneling of electrons interacting with vibra- 
tional modes. Theoretical studies of the coupling be- 
tween molecular vibrations and electronic states partici- 
pating in the tunneling have begun with the exact calcu- 
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lation of the single-electron transmission, in which the 
presence of the Fermi seas representing the leads has 
been essentially ignored.— i^i The single-electron trans- 
mission naturally exhibits resonances at energies corre- 
sponding to the vibration frequencies (these are often 
called "side-bands" ) . Such side-bands also appear in the 
local (on-molecule) single-particle density of states, com- 
puted in the presence of the leads j^^'^'^ There are claims in 
the literature^ii^ that these side-bands are reflected in, 
for example, the gate-potential dependence of the linear- 
response conductance. However, as has been emphasized 
by Mitra et ali^ and discussed in detail below, such side- 
bands cannot appear at zero temperature in the linear- 
response regime. We shall demonstrate their appearance, 
albeit weakly, at finite temperatures. 

Indeed, an electron crossing the molecular bridge may 
do so inelastically or elastically (with or without chang- 
ing the excitation state of the molecule). In the former 
case the electron will lose its phase coherence - a problem 
to which we will return in future work (see the discus- 
sion in Ref. '2^. Here, we concentrate on the structure of 
the conductance as a function of the bias voltage V and 
the gate potential, represented by the electrochemical po- 
tential ^ (when applying the latter is feasible). Clearly, 
at low temperatures only elastic processes and inelastic 
ones exciting the molecule are possible. The latter can 
happen only if the transmitted electron can supply the 
energy required for the molecular excitation. Focussing 
on a molecular vibration of frequency Wq, it is clear, then, 
that it can be excited only when the bias voltage V ex- 
ceeds hiOQ/e, namely, beyond the linear-response regime. 
Thiai^ will be confirmed by the detailed calculations be- 
low. 

The footprints of the inelastic processes appear in the 
differential conductance when plotted as a function of the 
bias voltage. (For analyses of the full counting statistics 
of a vibrating junction, see Refs. [l^.) This regime has 
been studied experimentally rather intensively. Theo- 
retically, it has been treated by employing a variety of 
methods and numerical techniques.— '^S'^i^'^ At low 
temperatures, the inelastic channel comes in when the 
bias voltage exceeds /lajg/e. This however does not nec- 
essarily imply an increase of the total conductance, since 
the elastic conduction channel might be modified as well. 
Indeed, interestingly enough, it has been observed that 
the "step" in the conductance at ^ = hujQ/e appears 
either as a decrease or an increase in the differential 
conductancei^iiiSiiSiiiS Theoretical work addressing this 
issue^i^ claimed that this behavior depends in a uni- 
versal manner on the bare elastic transparency of the 
jimction, T, such that the differential conductance steps 
upwards when T < 1/2, and downwards when T > 1/2. 
This claim has been refuted recently in a seminal paper 
by Egger and Gogolini^ We confirm their conclusion. 
Moreover, we find that the conductance steps downwards 
only in a narrow range of T, which becomes narrower as 
the ratio of the excitation energy Hloq to the bare reso- 
nance width Fr, increases. 



Another important aspect concerns the instabili- 
ties in the vibration modes possibly induced by the 
current In particular, Ref. [13 points out the in- 
applicability of the perturbation theory in the electron- 
vibration coupling once the noncquilibrium regime is 
reached. We show below that the step-like structure in 
the differential conductance at = fuoQ/e implies an- 
other type of breakdown of the perturbation theory. It 
turns out that the opening of the inelastic channel is 
inevitably accompanied by the appearance of logarith- 
mic singularities at the same bias voltage. Those are 
forced via the Kramcrs-Kronig relations and arc related 
to the singularities found by Engelsberg and Schrieffei— 
for bulk Einstein phonons. In this way we confirm the 
important (and seemingly, un-noticed) result of Mitra et 
al^ and Egger and Gogolin'^^: beside the step-like struc- 
ture, caused by the inelastic tunneling processes, the dif- 
ferential conductance develops a logarithmic singularity 
(at zero temperature, and to second order in the electron- 
vibration coupling) as the bias voltage crosses the vibra- 
tion energy. Near the threshold voltage V = Swp/e, that 
singularity dominates the differential conductance. 

It thus seems that there are several relevant issues 
in the theory of transport through a vibrating junction 
which are either still under debate or are not entirely 
clear. These concern the existence of side-bands, the 
dependence of the conductance on the junction trans- 
parency, the structure of the differential conductance 
near the opening of the inelastic channel, and the precise 
effect of the Franck-Condon factors on the resonances, in- 
cluding what happens when the resonance width exceeds 
the vibration frequency. Below, we give our answers to 
these questions and provide further physical interpreta- 
tions for them. In order not to obscure the basic physics 
by lengthy computations, we restrict ourselves to the sim- 
plest model, of a single resonance connected symmetri- 
cally to two leads and coupled linearly to a vibration. In 
addition, we apply lowest-order perturbation theory in 
the electron-vibration coupling. We believe that a com- 
plete analytical discussion of the outcome of this model 
will shed further light on the intriguing non-equilibrium 
behavior of the vibration-induced conductance. 

Section |TT] gives the Hamiltonian, and then expresses 
the current through the system in terms of the Green 
functions, which contain the contributions from the cou- 
pling to the vibrational mode. The detailed calculation 
of these Green functions is described in the Appendix. 
Section [Till presents the results for the conductance and 
for the density of states in the linear-response regime, 
while Sec. |IV] discusses the differential conductance at 
finite bias voltage (but zero temperature). Finally, we 
detail our conclusions in Sec. |Vl 



II. THE MODEL 

We consider the differential conductance of a small 
system, consisting of two leads connected together via a 
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"dot" . The two leads are assumed to be identical, except 
for being attached to reservoirs held at possibly different 
chemical potentials, /zl = /i + eV 12 and /ij^ = M ^ eV jl. 
When the electron is on the dot, it is coupled to a single 
harmonic oscillator of frequency loq . The Hamiltonian of 
this system is 



— ^Icad + ^dot + ^coup 



(1) 



The lead Hamiltonian is [using k(j>) for the left (right) 
lead, with the same lattice constant a = 1] 



with 



-2 Jcos k{p) 



(2) 



(3) 



The Hamiltonian of the dot is 



^dot = eocjco + f^o{b^b + i) + j{b + b^)clco , (4) 

where eq is the energy level on the dot, and 7 is the cou- 
pling energy of the electron (while it resides on the dot) 
with the oscillator. Below we often set h = 1. Finally, 
the coupling between the dot and the leads is described 
by 

Koup = J2 ^Mco + he) + ^ Vpiclco + he) , (5) 

k P 

with 



Hp) 



— JoSinfc(p) 



(6) 



(The wave functions on the leads are normalized assum- 
ing that each lead consists of N sites.) In Eqs. ([3]) and 
(ini), J is the overlap amplitude along the leads and Jq 
is the overlap amplitude between the leads and the dot 
(taken to be symmetric, for simplicity), all in units of 



energy. The operators cj, 



and c], 



and 



-p) cre- 



ate (annihilate) an electron on the dot, on the left lead, 
and on the right lead, respectively, while b"^ (6) creates 
(annihilates) an excitation of the harmonic oscillator, of 
frequency ujq. This model system has gained much the- 
oretical interest beforasee for example Refs. [H, [H, [l^, 

M [Ml, ii, [11 ii, il ml andii. 



details. The current entering the dot from the left lead, 



^T.^k{GU^)-G<,{u;)] 



X (-G<o(c.) + fL{co)[GU^) ~ GSo(cc')]) , (7) 



and that from the right one, Ird, is 



^RD — 6 



^E^p[^poM-Go<(c^) 
p 



duj 

«e / ^ro(^) 



X (-G<o(c.) + fR{Lo)[GU^) - G5o(^)] 



(8) 



In Eqs. O and ©, /l,k(w) = l/[e'^("-^i.,«) + 1] are the 
Fermi distributions in the two reservoirs, and Tq is the 
imaginary part of the self-energy Eq due to the coupling 
with the leads, Eqs. (|A11[) and (|A35|) . Obviously, current 
conservation requires /^^ + Ij^jj to vanish. Indeed, upon 
adding Eqs. (jT)) and ([5]) [and employing Eqs. (IA12P and 
(jA35p ] we find that current is conserved. Hence, the net 
current / can be obtained as the difference between the 
two currents, Ij^j^ and Ijid- This leads to 

I^e I ^r,{u;)[fr^{u;) - /«(c.)]ImGgo(u;) ■ (9) 



This well-known exact result in which the local density 
of states on the dot, ImGgg(a'), contains all of its dynam- 
ics, including the coupling to the oscillatoi\ is similar to 
the result as given, e.g., in Refs. [l^ 113, Hj, HI and[39l. 
It is customary to use the expression ([UJ also for non- 
linear transport. We remark that this is valid only for 
bias voltages that are not too large4S Equation © ne- 
glects the effects of the finite field on the system, such as 
the nonlinear screening, the induced changes in smd 
Jq, and the possibility, mentioned above, of "pumping" 
the molecule into higher states. Only when all the above 
finite-voltage corrections are neglected, does the Keldysh 
formulation justify using this result also in the nonlinear 
regime. 

The coupling with the harmonic oscillator affects the 
dot Green functions, Gqq and Ggp. In the absence of the 
coupling to the vibrations, the 'bare' dot Green function 
is given by 



A. The currents in the system 

The currents flowing in this system can be expressed 
in terms of the Keldysh Green functions: those on the 
dot are marked by the subscript 00, and the mixed ones 
are marked by the subscripts k(ji)0, see Appendix \X\ for 



^ - eg - Eg 



(10) 



Neglecting the frequency dependence of the self-energy 
due to the coupling with the leads, Eq, (this is the "wide- 
band approximation") and absorbing RcEq into eg, i.e., 
Eg — > e^jjg = + ReEg, the zeroth-order Green function. 
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Eq. ([T0|) . becomes 

^00 (^) 



(11) 



Expanding the Green function up to order 7^ yields 

g|;o = 4 + (4)'(A6o + sL) (12) 
= eo"o + (4)'fAi?±*ImSfJ , (13) 



where Ej^^ is the self-energy due to the coupling to the 
oscillator, Aeq is the shift of the energy, Eqs. (jAlsp and 

Aeo = -— / f^|e5o('^)Pro(^)[/L(^)+/;^(^)] , 
J 27r 

(14) 

and we have defined Ai? = Aeq + ReEg^. 

From the expansion Eq. it follows that the current 
can be written in the form— 



(15) 



where Iq is the current in the absence of the coupling 
with the oscillator, 

/d 
^r,{u;)[h{cj) - fj,iu;)]lmgSo{u;) , (16) 

/(.Q is the current involving the (real) shift in the resonant 
level (which depends on the frequency and the chemical 
potentials), 

X [{gS,iu^)f - {gUc^)f]AE , (17) 
and /jnco is the current involving the imaginary part of 

/d 
£ro(^)[/i(^)~/fl,(^)] 

X [{GSoi^))' + (e5o(^))']lniSL(^) • (18) 

Below we mainly consider zero temperature. (The ef- 
fects of a finite temperature on the linear-response con- 
ductance are considered in Sec. IIII CI ) Furthermore, we 
ignore the explicit dependence of Fq on w. At zero tem- 
perature, the zeroth-order current is 



Zn 



duj 



r n 



C.2+F2 



CFq / 

arctan— ^ — arctan— ^ 1 



2tt 



F / ' 
^ ^ 



(19) 



the current due to the effective shift in the resonance 
energy, Eq. (HZ]), is [see Eq. (fTT|) ] 



2 /-/J 






TT 



duj 



duj 



OJ 



(C.2+F2)2 



[AeoifJ-L^t^R) 



R.eI]go(w,/Xi,/i^)) , 



(20) 



and the current due to the imaginary part of the self- 
energy, Eq. (fTS]) . is 



eFn ^ uj^-n 
duj- 



° ImEg,(u;,/Xi,M^) . (21) 



27r ./„ "~(cc;2 + F2)2 



B. The zero temperature Green functions and 
self-energies 

The detailed calculations of the contributions to the 
Green functions due to the coupling with the oscillator 
are given in the Appendix. At zero temperature, Eq. 
(fT4l) becomes 



TTWn 



doj 



(w - e,,J2 + Yl 



(22) 



(It seems that this shift was overlooked in Ref. [S^-) In 
computing the explicit expressions of the currents and the 
conductances [see Sees. HIT] and [W], it is expedient^^ to 
measure the frequencies lo and w', as well as the chemical 
potentials Hl{b.) from e^.^^. We then find 



Aen 



7 



7 



■ E 

a=L.,R 



arctan 



(23) 



With the same notations, Eq. (jA36p gives 

imEhJc.) -7 (^ + ^^)2 + r2 

l-[/^(a>-^o) + /^(^-^o)]/2^ 



(^-^o)' + rg 



7'ro 



2^ I C/,, + , 



t^O - '^) , ©('^^ - Ma - ^o) 



(24) 



which reproduces the result of Ref. [S^- Clearly, 
ImE^Q(w) = unless w < /ii, — Wq and/or w > + i^o- 
Since this self-energy is required within an integral for 
which < UJ < jiL [see Eq. (pij) ]. its contribution to 
the current appears only when the bias voltage exceeds 
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TiLo^/e. Indeed, substituting Eq. (|24p in Eq. ((2T|) yields i^Q. The 8-function factor determines how much these 

Lorentzians contribute to the current. This reinforces the 
notion that I^^^^ is the current due to inelastic processes 
where a vibration quantum is given to or taken from the 
oscillator by the transmitted electron. As we discuss be- 
low, finite temperatures result in small contributions to 
the current /j^j-o even in the linear-response limit of zero 
bias voltage. 



duj 



1 



Ml 



duj 



(25) 



For Wq ^ Fq the integrand in Eq. (|25|) contains the 
two Forentzians shifted from the usual resonance by 



In a similar way, the real part of the self-energy is found 
from Eq. (jAMl) 



J 



+ 



2 



a; + 







-|- W| 







arctan 



a=L,R 



Fn 



To/ 1 

27r V (a; a;o)2 r2 



E In 

a=L,R 



n 



{UJ-II^+ Wo)2 (CJ - Wo)2 -I- F2 



E 1- 

a=L,R 



(26) 



again reproducing the result of Ref. Ha- [A simple interpretation of Eqs. ((24|) and ([26|) is given at the end of the 
Appendix, following Eq. (jA36|l ]. As mentioned, it is convenient to introduce [see Eq. (fT3|l ] the total energy shift 
which depends on the frequency and on the chemical potentials /i^ and ^r, 



AE{uj,iJ.L,f^R) = Aeo(/Xi,^;j) +Rei;;j„(w,^i,^fl) 

7^ro ^ln[(cj-^„-tJo)2/a;2] 



An 



a=L,R. 



{uj + uj„)^ + r2 



(27) 



where 



r 



[(:r-c.o)2+r2][(x + ..o)2 + F2; 



E In 

a=L,B. 



+ F^ - ^2) 

I 



ro)]z^ E ^r^tan^ 



„- 



a=L,R 



(28) 



The factor — 72/wQ, i.e., the polaron binding energy [see 
Eq. ([23]) ]. is independent of the frequency and of the 
chemical potentials. Therefore, we may safely absorb it 
in e^(,g and omit it from Eq. ((27|) . Inspection of Eq. ([27]) 
reveals that AE diverges logarithmically at w = /x^, ± . 
This divergence^ is dictated by the Kramers-Kronig re- 
lations once the imaginary part of the self-energy attains 
a discontinuity [see the discussion following Eq. ((24|) ]. 
The logarithmic divergence affects the conductance only 
in the nonlinear regime, and disappears in the linear- 
response one. However, the density of states is affected by 
these singularities even in the linear-response regime, see 
Sec. mil In any case, the logarithmic divergence implies 
that one should not ignore the frequency dependence of 
AE and absorb this energy in e^gg, as is sometimes done. 



III. THE LINEAR-RESPONSE REGIME 



A. Zero temperature conductance 



In the linear-response regime the bias voltage energy 
eV is the smallest energy, and at zero temperature the 
energy shift and the self-energy are required only a,t u) — 
fJ-L = IJ-R ^ where fi is the common Fermi energy 
of the leads (measured from the resonance energy e^.^^). 
Then, I-^^^ = 0, and I^^ [Eq. pp)) ] requires the energy 
shift AE{fj.,fi,^) = Afi{p, ^, fj.) [see Eq. ([77|) ]. which is 
a smooth function of fi. Thus, the only contribution to 
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the conductance from the couphng to the oscillator is 



2 CO 



temperature, the oscillator is in the ground state. Then, 
when the electron moves from, say, the left lead to the 



AE{fi,fi,^). (29) virtual state on the dot, the term 7(6 + 6'f)cJcQ in the dot 



In the linear-response regime the zeroth-order conduc- 
tance. Go, is 



2"g 



^ 



lin + r§ ' 

and we may combine Gq and Geo to obtain 
2tt 



-G 



lin [/i-Ai?(M,/X,/x)]'+r2 



(30) 



(31) 



Obviously, this expression is valid up to second-order in 
the coupling with the oscillator. 

AE/r,, 





(u„/r„=i/2 




<u„/r„=i 




FIG. 1: The function AE{fi, fi, fi), Eq. ^7}, for two repre- 
sentative values of the oscillator frequency, uq — CSFg (solid 
line), and uq = Fq (dashed line). Here 7 = Fq. 

The function AE(fi, fi, ^) is an odd function of /i, and 
its sign is opposite to that of /x, see Fig. [1] Therefore, 
in the linear-response regime, the effect of the coupling 
with the oscillator is just to reduce the conductance, 
except at resonance. This can be understood qualita- 
tively as due to the fact that the couplings to the two 
leads (Jo in our model) are renormalized downwards, 
to 0(7^), due to the same Franck-Condon-type factor. 
Therefore, the width of the resonance decreases, but its 
height, determined by the ratio of the two couplings, is 
unchanged.— The location of the resonance is not shifted, 
since fi—AE{^, fi, ^) = only at ^ = 0, i.e., —AE always 
moves away from the "bare" resonance energy. This be- 
havior is exemplified in Fig. [2l which shows the effect of 
coupling to the oscillator on the linear-response conduc- 
tance, for two values of the 'bare' width, Tq/ojq = 1/2 
(top) and 4 (bottom). 

Clearly, the relative narrowing of the resonance due to 
the vibrational mode decreases when the ratio Tq/ujq in- 
creases. Quantitatively, the 'renormalized' width of the 
resonance, F, is given by the solution of the equation 
r - AE{r,T,r) = Fq. As can be seen from Fig. [H 
AE{ii, n, n) becomes saturated at large /z; the solution 
for r increases towards Fq as Tq/luq increases. For a qual- 
itative understanding of this effect, we note that at zero 



Hamiltonian, Eq. (|4|), shifts the center of the oscillator 
motion by the order of 7/a;o (in units of the oscillator's 
zero-point displacement). However, this shift is fully real- 
ized only when the dwell-time of the electron on the dot, 
F,^"'^, is longer than the response-time of the oscillator, 
governed by ujq, i.e. Fq/wq ^ 1. In this limit, the cou- 
pling matrix element Jq will be reduced by the overlap 
integral between the shifted and the un-shifted oscilla- 
tor wave functions, which is of order cxp[— (7/wq)2/2], 
resulting in a relative narrowing of the resonancei^^ 

Indeed, at small 7/1^0 (^o that the perturbative ex- 
pansion is valid) and for small Fq/cjq we find that F/Fg 
approaches the Franck-Condon factor exp[— (7/0)0)^] 
(which becomes 1 — (j/ujfj)'^ in our order 7^ approxi- 



mation). This can be seen directly from Eqs. 
when ujq :§> {\^\, Fg}, Re'E^^{uj) is dominated by the 
first term in the large brackets on the second line, which 
is independent of the Fermi functions. Therefore, in 
this limit AE(fi, 11, fi) does not depend on the many- 
body effects contained in these Fermi functions, and 
the simple single particle Franck-Condon result is re- 
produced. Indeed, in this limit one has AE{ii, fi, ^) = 
Afi{fi,^,fi) fa —fi'j^/ujQ + 0(^Fo72/cjg), and therefore 
F « Fg(l — 72/u;o). In contrast, when Fg » ujq (bottom 
panel of Fig. the electron leaves the dot before the 
oscillator has responded to its presence, and the Franck- 
Condon blockade effect is much weakened. In our cal- 
culation, part of this blocking involves the Fermi func- 
tions on the leads [all the terms except the first in Eq. 
(PS)) ] . This dependence on the chemical potentials in the 
leads reflects the many-body effects on the leads, which 
seem to weaken the Franck-Condon blockade. So far, 
we have discussed the Franck-Condon narrowing only for 
zero temperature and zero bias voltage. However, the 
modified narrower shape of the resonances will also affect 
integrals over energy, causing apparently similar effects 
at finite temperatures and at a finite bias voltage. 

Another remarkable aspect is that at zero temperature 
the linear-response conductance exhibits no side-bands as 
a function of the gate voltage (modeled here by the com- 
mon /i measured from e^es)' when /i crosses the oscillator 
frequency. This has been emphasized in Ref. [TgI . con- 
trary to certain findings in the literature (see for exam- 
ple Refs. [2^ and |25[ ). Finite temperatures may generate 
small satellites, as discussed in Sec. IIII CI The absence 
of the side-bands in the linear-response conductance at 
the oscillator frequency, as the gate voltage is swept, may 
appear at first sight somewhat surprising. However, it is 
their appearance at zero bias voltage and zero temper- 
ature which is in fact un-physical. A structure in the 
linear-response conductance at ^ = i^o will mean that 
after passing, the electron leaves the dot in an excited 
state, even at zero temperature. As the electron begins 
and ends at almost the same energy, energy conservation 
does not allow it to excite the oscillator. 
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FIG. 2: The dimensionless linear-response conductance, Eq. 
pi[) . for two representative values of the ratio Fo/ci^oi 0.5 (top 
panel) and 4 (bottom panel). It can be seen that the conduc- 
tance (dashed line) is always smaller than that obtained in 
the absence of the coupling with the oscillator (depicted by 
the solid line).— 



B. The zero-temperature density of states 

The situation is very different when one looks at the 
local single-particle density of states on the dot, N{u!), 
given by 



-lmGoo(w) 



(32) 



This density of states is accessible, in principle, via lo- 
cal STM I — V measurements. We note that N{u}) is 
the nontrivial part of the integrand in the basic Eq. ^ 
for the current. Here we actually calculate it only at 
equilibrium, which is appropriate for the linear transport 
regime. For fiL — fJ-R = this quantity becomes 



1 



N{u;) = - 



ImS£o(^) 



TT lOJ - 



Ai?(u;,//,^)]2-f [ro+ImI]-Ja;)]2 



(33) 



where lmS^^{u}) is given by Eq. (l24)l . and AE{uj, fi, fi) is 
given by Eqs. ([26]), and Inspection of those 

expressions reveals that when fi — 0, i.e., the common 
chemical potential of the leads is aligned with the res- 
onance level on the dot, the density of states Eq. f33|) 
is even in the frequency, while at off resonance (where 
^ 7^ in our notations) it is not. In the first case, there 
will be equal weights for a hole (an electron) excitation 



corresponding to an excited oscillator and an electron (a 
hole). In the second, those weights are not equal. In par- 
ticular, when fJ, > 0, i.e., the common chemical potential 
of the leads is above the level on the dot, and there is 
more weight to the hole formation. 









wo/ro-1 








oj/Fo 



FIG. 3; The local density of states on the dot, Eq. 
as a function of the energy cu for /i = (top panel) and for 
jj, — 1.5Fo (bottom panel). The solid lines correspond to 
1^0 = ^0, while oJo = STo for the dashed lines (here , 7 = Fq). 
All the graphs should go to zero at lj = ±loq. 

The local density of states N{uj) is plotted in Fig. [3] 
at resonance, /x = (upper panel), and off-resonance, 
fl = I-SFq (lower panel). Both figures show structures 
around lu ~ /izbcjQ, which become smaller as luq increases. 
Each of these structures contains two ingredients: first, 
AE diverges logarithmically near uj = fi±ujQ, resulting in 
the vanishing of iV(a;) at these frequencies. Since these 
singularities are very narrow, the plots miss showing the 
actual vanishing of N{uj) at these points. The situation 
is somewhat more complicated for > and uj = fi — uJq 
(corresponding to the vicinity of wq/Fq = 0.5 for the full 
line in the lower panel of Fig. [3]). In that case, the energy 
difference ix' — Ai5(w, /i, fi) changes sign as one approaches 
the singular point, and therefore N{oj) first increases and 
only then decreases quickly to zero at w = /i — Wq. The 
plot picks up the initial increase, and misses the very 
narrow dip. 

The second effect arises from lmY,^^{uj), which modi- 
fies the width of the original resonance and creates the 
inelastic resonances. As can be seen from Eq. (|24p . 
this term contains contributions from two 'resonances', 



a,t LU — ±0Jq. However, the left (right) hand side reso- 
nance is included only for u < n — loq (uj > fj, + ujq) . For 
< luq, this causes a discontinuous increase in N{uj) 
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for uj below (above) n — oJq (// + Wq). The deep dips at 
Ltj = fi ± loq and the increased density of states beyond 
these energies create peaks in N{uj) at uj > ^ + ujq and 
at Lo < fj, — ujq, which can be identified as the side-bands 
(see top panel in Fig. [3]). For ^ > ujq, the behavior for 
Lu < ji — luq is more complex, but the general features re- 
main the same (lower panel in the figure). Note that our 
calculation shows only two such "side-bands" , since we 
work to second order in the coupling with the oscillator. 
Note also that the side-bands would not be as clear had 
we absorbed AE as a 'constant' in e,„,. 



C. Finite temperatures 

At finite temperatures, the linear-response conduc- 
tance is given by [see Eq. 
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FIG. 4: The linear-response conductance resulting from /jncoi 
for /SFq = 1, 3, 6 and 12 (increasing dash sizes). Here ujq = 
, and all energies are in units of Fg . 



duj(f3 



(e/5(c^-A0 + 1): 



, ImGgoH- (34) 



Using the expansion of Eq. (|13p , and the expressions in 
Eq. (|A36[) . and substituting /i/, — ^ji = /i for linear- 
response, yields 



and 



duj' 



UJ — UJ 



(coth^^l^ 



TT Uj'-^+Tl{u-Uj'Y 



^~^' tanh ^^^'~^'* 



(35) 



IV. THE ZERO TEMPERATURE 
DIFFERENTIAL CONDUCTANCE 



The differential conductance is the derivative of the 
current with respect to the bias voltage, V = {iij^—iip)/e 
[however, for finite bias voltage, note the discussion fol- 
lowing Eq. ©J. Differentiating Eq. with respect to 
V, the zeroth-order conductance is 



F2 
^ 



F2 
^ 



27r ^1 



(37) 



X ( coth 



stanh 



Pico 



SOJn 



M) 



(36) 



Figure [4] portrays the contribution of /jncoj Eq. (fT8| . 
to the linear-response conductance. We plot only this 
contribution, which arises from the inelastic processes, 
in order to exhibit the channel-opening due to the finite 
temperature. The contribution of I^^ is smooth, so it 
does not have drastic effects at finite temperature. Scal- 
ing /jj^^Q by /3Fq exp[— /Jwo], it is seen that the curves 
plotted for various temperatures approach an asymptotic 
limiting form for large /3Fg, exhibiting a reduction of the 
conductance near fi — and peaks slightly above (below) 
fi = luq (p = —ujo)- This structure could have been de- 
scribed as having side-bands; however, the peaks decay 
exponentially (as exp[— /Jwq]) a-t low temperatures. This 
factor arises directly from the low temperature behavior 
of the large brackets in Eq. , and is also understand- 
able intuitively: the side-bands can contribute only if ex- 
citations by the oscillator energy Hujq are allowed. Those 
appear with the Boltzmann factor exp[— /^wq]. As the 
temperature increases, the structure portrayed in Fig. [4] 
broadens and gradually becomes smeared. 



Similarly, differentiating Eq. (|^D|) gives 



where 



9 ^CO 9 ^ CO ' 9 ^CO ' 



(38) 



9 ^ CO 



and 
27r 



E 7^^x^Ai?(^ = /i„,/ii,/i«) , (39) 



a=L,R. 



2 '-'co — 



duj- 



dAE{Lj) 



Using Eq. (|27p . we have 

dAEicu) _-/^r„f 1 {lu - Mfl) Vc^o 



(40) 



In \ul + 



A rather lengthy computation yields 



(41) 
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27r 



+ 



+ 



E 

a=L.R 



2r3 r 





(A.i + r2)(^|^ + r§) 



47r 



47r 
47r 



^ 



/;2 -I- r2 



To -/^lW +^0) 



In 



(42) 



where A/z is given by Eq. (|28|) . and we have defined 
1 /(x + y)(y2+r2) 



J/ In 







r2 + r2 

-^^ arctan — 



(43) 



Also, a marks the lead which is not a. One observes that 
when fi]^ = fiji, then G^^ is fully given by only the first 
sum in Eq. (|42p . reducing to the linear-response result 

Finally, the differential conductance resulting from the 
current li^co^ Eq. ([25]) . is 



-G = 

^inco 



^e(^, 



^ 



(44) 



At the threshold bias voltage, eV = IJ'j^ — M/j — Wq, this 
contribution jumps from zero to 



27r 



^ 



(A*i + rg)2(^?,, + r2)2- 



(45) 



Since at eV — ujq the chemical potentials are fj.^ = fj, + 
uj()/2 and /i^ = fi — uJo/2, it follows that the conductance 
jumps downwards when the common chemical potential 
of the leads (measured from the resonance level) is in the 
range 



max[0, {ujj2f - Vl] < < {ijj2f + E^ 



(46) 



Note that the range in which the conductance jumps 
downwards is shrinking as loq/Tq becomes larger. Since 
the "bare" elastic transparency of the junction is given 
by T = Eg/(^2 + Eg), the condition ([^ can be put in 
the form 



1 



2 + (C.0/2E0) 



2 = ^1 < ^ < ^2 = min[l, {2rjuJof] 



(47) 



One notes that the lower border-line transparency^^^ 71 
reaches the "universal" value 1/2 (with T2 = I) only^ in 
the limit of a very broad resonance, Eq ^ Wg, where the 
effects of the vibrational excitations are smeared within 
the original resonance. As the ratio cjq/Eq increases, both 
7i and T2 decrease and approach each other, so that the 
region with a negative step in Gj^j.^ narrows down in this 
physically relevant region. 



Eor V > 0, the logarithmic divergence in G^^ arises 
only due to the last term in Eq. (|42p . Near eV = Aiwq, the 
coefficient of this term contains the factor Eg — MlM/? — 
Eg -I- (a;g/2)^ - which is positive for all T > Ti. Since 
the argument of the log is small near eV = Itioq, this 
implies a negative divergence of this term in this range 
(and a positive one for T < %). Interestingly, the log- 
arithmic term does not change sign at T2, although the 
step in Gijjj.Q does change sign there. We note that the 
Kramers- Kronig relation, relating these two singularities, 
applies only to the real and imaginary parts of the self- 
energy, and not to the corresponding contributions to the 
differential conductance. 



Figure [5] shows the total conductance G = Gq + G^^ + 
Gijjj.Q, as well as the two separate contributions from the 
coupling to the oscillator G^^ and G-^^^^^ for Wg = 3Eg and 
for five values of T. The plot for G^^ does not contain 
the first term in Eq. (UH), which was incorporated into 
Eq. ([57)1 by the replacements /z^ n^ — A^{ijl^, Hj^, Hfj). 
Clearly, there are no visible singularities when T = 7i, 
where the coefficients of both the logarithmic term and 
the discontinuity vanish (there remain effects for higher 
derivatives of the current). The former singularity sur- 
vives at T = 72, where the discontinuity vanishes (al- 
though more steeply than near 7i). The logarithmic di- 
vergence is indeed positive for T < 7i. Also, the mag- 
nitudes of both the jump and the logarithmic divergence 
are large at large bare transparencies T, and decrease 
with decreasing T. It should be kept in mind that the 
apparent divergence in G^^ results from our expansion in 
powers of 7, which breaks down very close to the thresh- 
old V — hw^/e. 
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FIG. 5: The total differential conductance (left panel) and the contributions to the conductance from the coupling to the 
vibrational mode, Geo (middle panel, containing the logarithmic divergence) and G^^^^ (right panel, containing the discontinuity; 
Ginco = for eV < huiQ, and the jumps in the various conductances at the thresholds are seen as the lapses in the curves) 
as function of the bias voltage, for ujq = SFg, and five different values of the zeroth-order transmission T — 0.2, Ti = 
4/17, 0.3, 72 = 4/9, 0.9 (dash sizes increase with T, appearing in increasing order at small V in the left panel). Here 7 — Fq. 



V. SUMMARY 

Although several of our formal results have already 
been obtained by Mitra et al}^ and by Egger and 
Gogolin,-^ our paper has used these results for a criti- 
cal discussion of several physically relevant issues which 
have been debated in the literature. Specifically, we have 
obtained and highlighted the following points: 

(A) At zero temperature, the resonance peak in the 
linear-response conductance always narrows down due 
to the coupling to the vibrational mode. However, this 
narrowing down is given by the Franck-Condon factor 
only for narrow resonances, and when one may ignore 
the Fermi statistics of the electrons on the leads. When 
Fq ^ oJq, the electron dwell-time on the dot is short and 
therefore the relative narrowing is much smaller. 

(B) Contrary to claims in the literature, the linear- 
response conductance does not exhibit any side-bands at 
zero temperature. Small satellites, of order exp[— /Swg], 
do arise at finite temperatures, where the excitation of 
the vibrational mode becomes possible. 

(C) The coupHng to the vibrational modes does show up 
in the single-particle density of states, which exhibits two 
singularities at the frequencies uj — fj, ± loq which corre- 
spond to the opening of the inelastic channel in which 
the vibrational mode remains excited. (We find only two 
singular points, because we expand the results only up 
to second order in the coupling to the vibrational mode.) 
These include discontinuities, due to the imaginary part 
of the self-energy, and logarithmic singularities, due to 
the real part of the self-energy. The latter result in deep 
dips in the density of states around each threshold, cre- 
ating apparent side-bands at frequencies which exceed 
these thresholds. Although a logarithmic singularity im- 
plies the inapplicability of the perturbative expansion 
very close to the inelastic thresholds, the predictions of 
dips and satellites in the density of states can probably 
be trusted out of these narrow regions. 

(D) The same singularities also generate discontinuities 
and logarithmic divergences in the differential conduc- 



tance at and around the thresholds eV — zLHujq. The 
signs of the discontinuities are usually positive, but they 
become negative within a finite range of the bare elastic 
transparency of the junction, shrinking progressively as 
Wq/Fq is increased. The "universal ratio"^^^ is obtained 
only-^^ in the limit ojq /Fq <C 1 , in which the state with an 
electron on the dot and the vibrational mode excited is 
not well-defined. In contrast, the logarithmic divergences 
remain negative for a rather broad range of bare trans- 
parencies, indicating the breakdown of our perturbative 
expansion very close to the inelastic thresholds. 
(E) Contrary to some claims in the literature, our results 
are quite different from those based on the single-electron 
transmission, which ignore the Fermi seas in the leads. 

It would be useful to test some of these predictions in 
experiments or in other theoretical calculations. 
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APPENDIX A: THE KELDYSH GREEN 
FUNCTIONS 

We obtain the Green functions of our system by solving 
the Dyson equations up to second order in the coupling 
7. In this procedure, we use the following relation^ for 
the lesser product of two Green functions, 

{AB)< ^A''B< +A<B'' , (Al) 
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and similarly for the greater product of two Green func- 
tions, denoted by the superscript >. Here, 



(A2) 



Note that when the operators a and b are identical, 
and are purely imaginary. Another property of these 
Green functions (for general a and 6) is 



G<H=z J dte^-\b^a{t)) , 
G> (c) = -^ / dte^-'iaitp) 



(A3) 



where G*" (G") is the retarded (advanced) Green func- 
tion, 



Jo L J + 

cO 



I I e 

' — oo 



i(w-iO+)t , 



a(t),6t ) 



(A4) 



For brevity, the frequency w does not appear explicitly 
in most of the equations below. 

The Dyson equation of the Green function on the dot. 
Goo, reads 



Goo = .90 (l + VkGko + ^p'^po + ^'^1 



QQO 



(A5) 



Here, go is the free Green function of the dot (in the 
absence of the coupling with the harmonic oscillator and 
with the leads), i.e., go — — eo)^^. The other Green 
functions in Eq. (|A5[) are those mixing the leads and the 
dot operators, 

Gk(p)o =^ ^k{p)''4) ^ > (A6) 
and the one mixing the dot and the oscillator operators. 



G, 



OQO 



(A7) 



In the notations <C; the first (second) operator (or a 
product of operators) is the operator denoted by a (b) in 
Eqs. (1X2]) and 

The Dyson equations of the Green functions (|A6|) are 



Gk{p)o — 9k{p)^k{p)Goo 



(A8) 



Here gk{p) is the free Green function of the left (right) 
lead, 



1 



^ - %P) + ''0 



5fe(p) 



9k{p) 

9k{p) = (fffc(p) - 9k{p))fL{R){^) 



2m5{e 



Hp) 



(A9) 



where fL{R) is the Fermi distribution of the left (right) 
reservoir. As mentioned above, we assume that the two 



leads are identical except for their different Fermi func- 
tions. It therefore follows that 



^fcGfeo + Z ^pGpo — SqG, 



00 ) 



(AlO) 



where Eq is the self-energy due to the coupling of the dot 
with the leads, 



, - e, + iO+ ^ w - e„ + iO+ 

J, K p f 



^ in — 



and 



^< _^o ^0 



(/i^ + /l) 



(All) 



(A12) 



Thus, the Dyson equation (|A5|) of the dot Green function 
becomes 

(go-i-E;;)G;;o-i + 7G5Qo > (A13) 



and 



(go-i-ES)G<o = E<GSo+7G, 



< 

OQO • 



(AM) 



In particular, the dot Green functions in the absence of 
the coupling with the harmonic oscillator, Qooi are 



^00 ~ [ LO — En — T, 



^0 J ' 



»oo — yoo'^o »oo — ; 



(% - G'oo) ■ (A15) 



The self-energy coming from the coupling with the har- 
monic oscillator results from the Green function Goqo, 
Eq. (|A5|) . Its Dyson equation reads 



G, 



OQO 



{b + b^) 



E ^fcGoQfe + E ^pGoQp 



(A16) 



where 



and 



G, 



OQfc(p) 



GoQOQ =« Co(^' + b^y, {b + b^H » . (A18) 
It is straightforward to obtain 

E ^kGoQk + E ^pGoQp = SqGoqo j (A19) 
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and 



Thus we find from Eq. (jAl6|l tliat 

(5o-^-s;s)G;5Qo = (fc + ^'^) + 7GSQOQ 



and 



OQOQ 



(A20) 



(A21) 



(A22) 



Inserting the result (|A2ip into Eq. (|A13p gives that 
the retarded Green function on the dot, up to second 
order in the couphng 7, is 



where we have defined 



OQOQ • 



(A23) 



(A24) 



An analogous result holds for the advanced Green func- 
tion. Using this result and Eq. (|A22|) in Eq. (|A14|) 
yields the Keldysh Green function on the dot (again, up 
to second order in 7), 



^ho )'-^oo ' 



with 



OQOQ • 



(A25) 



(A26) 



It is hence found that the coupling with the harmonic 
oscillator modifies the dot Green function in two ways. 



J 



Firstly, it adds the term Sj^^ = j'^Gqqoq to the self- 
energy. This contribution is calculated below. Secondly, 
it shifts the resonance level by the amount 



2172 



(A27) 



This result is found by employing perturbation theory. 
To first order, the oscillator wave functions can be written 
in the form 



^ n' 



n' 



t (n\b + b''\n') , . 



and consequently the diagonal (n'n') matrix element of 
h + h^ is 



27 



{n\b\n'){n'\h'^\n) + {n\h'^\n'){n'\b\n) 



27 



-Cr,C| 



o'-o 



(A29) 



The average of cJcq is needed to zcroth order in the cou- 
pling with the oscillator, and therefore is expressed in 
terms of the Green function (jAlsp . leading to Eq. (IA27p . 

It remains to compute the Green function Gqqoq , Eq. 
(jAlSp . As we work up to second order in the coupling 
7, it is enough to find this function in the absence of the 
coupling to the oscillator. At this order, the electron op- 
erators and the oscillator operators are decoupled. For 
example, using the definitions of the Keldysh Green func- 
tions, Eqs. ((J2|l and (|X4| . 



GoqoqH = I dte^-\{m+b\mb + b^)){c,{t)cl) + {{b + b^Mt) + bHt))){clc,{t)) 

pOO nOO pOC nOO 

Jo Jo Jo Jo 

Therefore, upon carrying out the time-integrations, we obtain 

GoqoqM-*; ■^1^oo(^)(^_ , 



G^oi^') 



' ±iO+ lu + lUq-lu' ± iO+ 



! + UJq — LU' i: iO+ to — LUq 



1 )] 

y ±M+J i 



(A30) 



(A31) 
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Similarly, 



Go<qoqM=* J dte^-'{{b + b^){b{t) + bHm{clc,{t)) 
= t{bb^) f dte^(-+-o)*}(ctco(0)+z(fetfo) / o)*)(4co(t)) , 



(A32) 



and consequently 

(A33) 

where Qqq is the Green function on the dot in the absence 
of the coupling with the oscillator, see Eqs. (|A15|1 . It is 
easy to check that at zero temperature and at equilibrium, 
Eq. (jA32p reduces to the usual diagrammatic expression, 
see e.g., Ref. H 

In order to present explicit expressions for the self- 
energy due to the harmonic oscillator, we use [see Eqs. 



(IIl2| and jUSI] 

^o>oM = *ieSoHProMf/LM + /flM-2) , (A34) 



where we have denoted 



EgH-s]5(c.) 

2i 



It follows that [see Eqs. (|X24l) and (IMTj) ] 



(A35) 



rli )' 



UJ + cUq — to' ± iO+ 



(A36) 



It is instructive to interpret Eq. (jA36p in the simple 
equilibrium case where /l(w') = /fl(w') = /(w') as the 
change, within second-order perturbation theory in 7, 
of the energy of an electronic state at energy w, due 
to all other states, at a running energy u'. The term 
|^oo('^')PrQ((jj') appearing before the large brackets is 
just the density of the latter states at zero-order in 7. 
For the imaginary part of Tt^^^\uj), the first term in the 
large brackets is due to real transitions occurring by ex- 
citing the oscillator (intensity proportional to {bb'')) and 
going to ui' = (jJ — (jJq with the blocking factor 1 — /(w'), or 



by absorbing a 'phonon' (intensity proportional to (b'^b)) 
and going to uj from the same w', now with an initial 
population f{uj'). The real part of Yilj-^\uj), given by 
the principal part of the integrals, is just the correspond- 
ing perturbation-theory energy shift. Obviously, these 
real and imaginary parts satisfy the Kramers-Kronig re- 
lationships. The second term in the large brackets is 
likewise understood as involving transitions to the state 
uj' = uj + ujq (for the imaginary part) or to the states 
around it (for the real part). 
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